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ABSTRACT 


Anti-mine  blast  deflectors  are  now  used  as  a  standard  countermeasure  to  re¬ 
duce  landmine  effects  on  a  vehicle’s  occupants,  increasing  their  chances  of  survival. 
The  blast  deflectors  are  submitted  to  an  impulse  from  the  soil  ejecta  and  detona¬ 
tion  products  resulting  from  a  mine  explosion.  This  memorandum  presents  algebraic 
equations  to  estimate  the  total  impulse  imparted  to  these  deflectors.  The  proposed 
equations  are  suitable  for  use  in  a  spreadsheet  program  and  can  be  used  by  engineers 
and  scientists  to  develop  anti-mine  protection  systems. 


RESUME 

Des  deflecteurs  de  souffle  antimine  sont  maintenant  utilises  de  fagon  stan¬ 
dard  pour  diminuer  les  effets  de  mines  sur  les  occupants  d’un  vchicule,  augmentant 
ainsi  lcurs  chances  de  survie.  Les  deflecteurs  de  souffle  sont  soumis  a  l’impulsion 
causee  par  lc  sol  et  aux  produits  de  detonation  resultant  de  l’explosion  d’une  mine. 
Ce  memorandum  presente  des  equations  algebriques  pour  estimer  rimpulsion  totale 
transmise  a  ces  deflecteurs.  Ces  equations  peuvent  etre  programmees  dans  un  chiffrier 
electronique  et  utilisees  par  les  ingenieurs  et  scientihques  qui  mettent  au  point  les 
systemes  de  protection  antimine. 
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EXECUTIVE  SUMMARY 


Landmines,  also  called  anti-tank  mines,  vehicle  mines  or  simply  mines,  have 
been  encountered  in  recent  peacekeeping  missions  in  Cambodia,  Somalia  and  Bosnia. 
Even  though  the  Canadian  Forces  have  acquired  some  mine-proof  vehicles  from  South 
Africa,  there  is  a  need  to  develop  anti-mine  protection  systems  for  the  Canadian  B- 
fleet,  ranging  from  the  light  support  vehicle  wheeled  (LSVW)  to  the  heavy  engineer 
support  vehicle  (HESV).  To  help  meet  this  requirement,  the  Defence  Research  Es¬ 
tablishment  Valcartier  has  undertaken  a  defence  R&D  project  to  better  protect  the 
occupants  of  soft-skinned  vehicles  from  landmine  effects. 

The  capability  of  estimating  the  impulse  transmitted  to  a  vehicle  is  essential 
to  predict  or  simulate  the  motion  of  the  vehicle  and  its  occupants  subjected  to  a 
landmine  explosion.  This  memorandum  provides  algebraic  equations  to  estimate 
the  total  impulse  on  vehicle  plates  acting  as  blast  deflectors — an  important  counter¬ 
measure  to  deflect  mine  effects  away  from  the  vehicle’s  occupants.  These  equations  are 
applied  in  a  series  of  examples  to  obtain  an  estimate  of  the  magnitude  of  the  expected 
impulses  and  to  study  the  effects  of  some  changes  in  the  mine/target  geometry  on 
the  resulting  impulses. 

The  impulse  produced  from  a  landmine  explosion  is  mainly  the  result  of  deto¬ 
nation  products  and  soil  ejecta  impacting  the  target  at  high  velocity.  The  U.S.  Army 
Tank-Automotive  &  Armaments  Command  developed  in  1985  an  empirical  equation 
to  predict  the  vertical  component  of  the  specific  impulse  (impulse  per  unit  area)  for 
buried  charges  with  a  pancake  shape.  That  vertical  component  of  the  specific  im¬ 
pulse  is  used  in  this  memorandum  to  derive  algebraic  equations  for  the  total  impulse 
on  horizontal,  oblique  and  vertical  blast  deflectors.  An  expression  for  the  maximum 
impulse  to  be  expected  for  a  given  standoff  distance  is  also  presented. 

These  equations  will  be  used  in  the  development  of  efficient  anti-mine  protec¬ 
tion  systems  at  DREV.  However,  more  research  is  required  to  combine  the  impulse 
on  blast  deflectors  to  the  impulse  on  a  wheel  for  an  explosion  under  a  tire. 
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1  INTRODUCTION 


Since  the  last  World  War  and  the  creation  of  the  United  Nations,  Canada 
has  been  involved,  to  various  degrees,  in  all  UN  peacekeeping  missions.  Landmines, 
also  called  anti-tank  mines,  vehicle  mines  or  simply  mines,  have  been  encountered 
in  recent  peacekeeping  missions  in  Cambodia,  Somalia  and  Bosnia.  While  most  UN 
missions  do  not  put  the  lives  of  the  peacekeepers  at  immediate  stake,  these  three 
cited  missions  did  indeed  put  soldiers  in  hazardous  situations  and  good  anti-mine 
protection  systems  should  have  been  required  for  most  of  the  vehicles. 

In  Bosnia,  53%  of  landmine  accidents  occurred  while  the  vehicles  hit  mines  on 
“cleared  routes”  (see  [1]),  showing  that  in  a  mine-infested  region,  one  must  always 
be  protected  from  landmine  effects.  Even  though  the  Canadian  Forces  have  acquired 
some  mine-proof  vehicles  from  South  Africa,  there  is  a  need  to  develop  anti-mine 
protection  systems  for  the  Canadian  B-fleet,  ranging  from  the  light  support  vehicle 
wheeled  (LSVW)  to  the  heavy  engineer  support  vehicle  (HESV).  To  help  meet  this 
requirement,  the  Defence  Research  Establishment  Valcartier  has  undertaken  a  defence 
R&D  project  to  better  protect  the  occupants  of  soft-skinned  vehicles  from  landmine 
effects  (see  [2,  3]). 

To  predict  or  simulate  the  motion  of  a  vehicle  and  its  occupants  submitted 
to  a  landmine  explosion,  one  needs  the  capability  of  predicting  the  impulse  trans¬ 
mitted  to  the  vehicle.  This  technical  memorandum  provides  algebraic  equations  to 
estimate  the  total  impulse  on  vehicle  plates  acting  as  blast  deflectors — an  important 
counter-measure  to  redirect  the  mine  effects  away  from  the  vehicle’s  occupants.  These 
equations  are  applied  in  a  series  of  seven  examples  to  obtain  an  estimate  of  the  mag¬ 
nitude  of  the  expected  impulses  and  to  study  the  effects  of  various  changes  in  the 
mine/target  geometry  on  the  resulting  impulses. 


These  algebraic  equations  will  be  helpful  to  the  engineers  who  need  to  estimate 
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the  loads  for  the  design  of  blast  deflectors,  and  to  the  scientists  who  develop  anti-mine 
protection  systems. 

To  describe  the  mine/target  geometry,  a  fixed  system  of  coordinates  (the  ob¬ 
server’s  view)  is  first  defined  in  Chapter  2. 

Westin  et  al.  in  [4]  have  developed  an  empirical  equation  for  the  specific  im¬ 
pulse  transmitted  to  plates.  Chapter  3  proposes  an  approximation  of  that  empirical 
equation  from  which  derived  algebraic  equations  for  the  total  impulse  are  found  and 
presented  in  Chapter  4  for  horizontal  deflectors;  in  Chapter  5  for  oblique  deflectors; 
and  in  Chapter  6  for  vertical  deflectors. 

This  work  was  carried  out  at  DREV  under  project  work  unit  2bal5,  Surviv¬ 
ability  and  Lethality  Assessment  and  Modeling  Software ,  between  January  and  March 
1998. 


2  REFERENCE  SYSTEM 


The  origin  O  of  the  xyz-axes  is  located  at  the  center  of  the  mine,  buried  at 
a  depth  5,  as  shown  in  Fig.  1.  The  horizontal  plane  is  the  xy-plane  and  the  z-axis 
points  towards  the  sky. 


3  FUNDAMENTAL  MODEL 


The  impulse  produced  from  a  landmine  explosion  is  mainly  the  result  of  det¬ 
onation  products  and  soil  ejecta  impacting  the  target  at  high  velocity  (see  [4]).  The 
loading  generated  from  the  impulse  is  characterized  by  a  high  pressure  during  a  short 
period  of  time.  The  duration  is  much  shorter  than  the  natural  period  of  the  vehicle 
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FIGURE  1  -  The  Coordinates  System 

structure,  and  so  the  exact  shape  of  the  pressure-time  loading  is  not  as  important  for 
the  calculation  of  the  vehicle’s  response  as  its  integral  over  time.  The  integral  of  the 
pressure  loading  over  time  is  called  the  specific  impulse. 

In  [4],  Westin  et  al.  proposed  an  empirical  equation  to  predict  the  vertical 
component  of  the  specific  impulse  (impulse  per  unit  area)  for  buried  charges  with  a 
pancake  shape.  The  vertical  component  of  the  specific  impulse  is  fundamental  for  the 
characterization  of  the  total  impulse  since,  as  will  be  shown  in  Chapters  5  and  6,  the 
other  components  are  easily  derived  from  the  vertical  component. 

3.1  Vertical  Component  of  the  Specific  Impulse 

Let  P  =  P(x,  y,  z)  be  a  point  on  a  horizontal  blast  deflector,  as  shown  in 
Fig.  2.  The  empirical  expression  for  the  vertical  specific  impulse  iv  at  P  is,  from  [4], 


[Pa-s]  (1) 
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where 


and 


z  =  standoff  distance  of  P  to  the  center  of  the  mine  [m] 
d  =  \Jx2  +  y2  =  lateral  distance  to  center  of  the  mine  [m] 
6  =  burial  depth  to  center  of  the  mine  [m] 


p  =  soil  density  [kg/ 


nr 


E  =  energy  release  in  explosive  charge  [J] 
A  =  cross-sectional  area  of  the  mine  [m2l 


c  = 


-5/4  ^3/s  tanh((2.2f) 


3/2 


m 


-ii 


( - 

1 

S 

- v  1 

Mine  (E,A)  f 

Blast 
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FIGURE  2  -  Horizontal  Blast  Deflector 

Because  Eq.  (1)  is  an  empirical  one,  it  should  be  used  only  when  those  four 
criteria  are  all  met: 


0.106  <  -  <  1.00 

z 


(2) 


E/A 

6.35  <  -A—  A  150.0 


pc2  z 


(3) 


5 


0.154  < 


Va 


<  4.48 


(4) 


0  ^  ^  19.3  (5) 

The  parameter  c  in  Eq.  (3)  is  the  seismic  P-wave  velocity  in  the  soil  expressed  in 
[m/s]. 


The  function 


u/-  ,n  def  ( tanh(0.9589  (  d) 

c  d 


3.25 


in  Eq.  (1)  is  not  defined  at  (  d  =  0;  when  (  d  — >  0,  /(£  d)  goes  to 


(6) 


limo/(C  d)  =  (0.9589)3'25  =  0.8725. 


Then  for  small  values  of  Q  d  {e.g.  (d  <  0.05  )  Eq.  (1)  for  the  vertical  component  of 
the  impulse  becomes 


iv(x,  y )  —  0.1180 


(7) 


3.2  Modified  Empirical  Function 


The  total  impulse  on  a  blast  deflector  is  given  by  the  integral  of  the  specific 
impulse  over  the  total  area  of  that  blast  deflector. 


=  //  iv(x,  y)  dx  dy 


[N-s  . 


An  algebraic  equation  for  the  total  impulse  is  sought  in  this  memorandum.  To 
integrate  analytically  Eq.  (1)  is  made  impossible  by  the  form  of  the  function  f((d) 
defined  with  Eq.  (6).  The  solution  proposed  in  this  document  is  to  approximate  /(/  d) 
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by  a  function  for  which  an  analytical  solution  of  the  integral  of  the  specific  impulse 
over  the  entire  plate  surface  exists.  The  function  g{Q  d)  defined  as 

/( C  d)  «  g((  d)  =f  K!  (l  +  k2  (C  df )  e"K3  {Cd)\  g{ 0)  =  lim  f((  d)  (8) 

TM 

meets  this  requirement  and  its  parameters  were  sought  using  the  Mathematica 
function  NonlinearFit  (see  [5]).  The  best  values  for  the  fit  are 

ki  =  0.8725,  k2  =  0.04837,  k3  =  0.8917.  (9) 


Fig.  3  compares  the  functions  /  and  g  and  the  absolute  error  f  —  g  is  plotted  in  Fig.  4. 


FIGURE  3  -  The  Approximation  of  /(C  d)  by  g((  d) 


0  1  2  3  4  5  6 

;d 


FIGURE  4  -  The  Absolute  Error  f((d)  -  g((d) 
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4  HORIZONTAL  BLAST  DEFLECTORS 


4.1  Total  Impulse 


The  total  impulse  iv  on  a  horizontal  blast  deflector  is  given  by  the  integral  of 
the  specific  impulse  iv(x,  y )  over  the  entire  plate  surface.  If  the  projection  of  the  blast 
deflector  on  the  xy-plane  is  as  seen  in  Fig.  5,  then  the  total  vertical  impulse  on  the 
deflector  is 

rx i  ryi 

W  =  /  /  iv(x,y )  dydx 

Jx0  Jyo 

=  0  1352  (1  +  ^) 


o  E  rx i  rv i 

\  -  /  /  f  {(d)  dydx.  (10) 

V  Z  Jx0  Jy0 


FIGURE  5  -  Horizontal  Blast  Deflector  over  the  xy-Plane 


4.2  Approximate  Solution 


There  are  no  known  analytical  solution  for  the  integral  in  Eq.  (10).  This 
integral  can  be  numerically  integrated  but  the  process  is  time  consuming  and  does 
not  generate  an  algebraic  equation  for  the  total  impulse.  The  solution  proposed  is  to 
use  the  approximation  (8)  to  replace  f((d)  by  g((d)  in  Eq.  (10). 


The  vertical  impulse  on  the  deflector  can  now  be  rewritten  as 


where 


iv  -  *0  (.S.  +  s2)  (l  +  \j^f 


k0  =  0. 1352  =  0.1180 


Si 


XI  ryi 


,-«3  (C^)2 


dy  dx 


xo  Jy o 


7T  AX  Ay 

4^3  C2 


*5*2 


(C  d )6  e  K3«d)2  dy  dx  =  16^4C2  E  a i 
Ax  =  erf(y^xiC)  -  erf(yTc3  xo  C) 

Ay  =  erf (v^ 2/i  C)  -  erf(v^2/o  C) 


fX  1 

/lo 


yi 

yo 


(11) 

(12) 

(13) 

(14) 

(15) 

(16) 


-4i  =  ySiiiCe-'‘3';<a 


X 


12  v/^3  C  ( e~K3y ^ 


_e~K3  yg  <2 


(3  2/1  +  k3  x2  2/1  C2  +  k3  y\  C2) 

2  y-2  ,  3  a 2\\ 

2)1  1/oC  +«3//0C  JJ 

v^vr  Ay  (33  +  16  k3  x\  £2  +  4  ^  x4  C4) 


(IT) 


A2  =  3  yTc3  x0  C  e_K3  x°  ^  x 

4  C  (3  l/o  e“K3  ^  c2  -  3  2/1  e“K3 ^  +  k3  1 /03  C2  e""3  *°2  <2 

—^3  2/?  C2  e*32^2)  +  11  v^Ay 


(18) 


A3  =  -4v^4xE5AJ/e-K33:oC2 


(19) 


A4  =  4  /c|  x3  C3  e"K3  *0 f2  [3  C  (2/0  e"K3  ^  -  2/1  e“K3  e)  +  4  ^  A4 


(20) 
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A5  =  v^Aa 


33  y/Ki  y0  C  33  yfni  y\  (  16  /t32  y„  C3 


=  2/q  f2 


3«;3  y\  C2 


+ 


3*3  2/q  C2 


IQnjylC  4k32^C5  4/cf^C5 


=«3  y  {  C 


+ 


=  «3  C 2 


=  K3  3/i  C2 


+  24  \/7r  Aj 


(21) 


The  function  erf(-)  in  Eqs  (15)  and  (16)  is  called  the  error  function.  The 
appendix  gives  the  definition  of  the  error  function  and  presents  an  algorithm  for  its 
evaluation. 


Example  1  To  compare  a  numerical  solution  of  Eq.  (10)  to  the  proposed  algebraic  solution  (11),  consider  the 
explosion  of  a  TMA-3  anti-tank  mine  under  a  vehicle  equipped  with  a  horizontal  blast  deflector  installed  40  cm  above 
the  ground.  The  deflector  area  is  1  m  by  1  m  and  is  centered  above  the  mine.  The  TMA-3  mine  is  26.5  cm  in  diameter 
and  8  cm  in  height,  contains  6.5  kg  of  cast  TNT  and  is  buried  under  3  cm  of  soil  of  density  1.6  g/cm3.  Then 


s  =  0.40  +  0.03  +  ^  =  0.47  [m] 

•To  I/O  •-  -  2  [m] 

XI  =  yi  =  \  [m] 

5  =  0.03  +  2^8  =  0  07  [m] 
p  =  1600  [kg/m3] 

E  =  6.5  [kg]  x  4.516  [MJ/kg]  =  2.935  107  [J] 


A  = 
C  = 


7r  X  0.2652 


=  0.055  [m2] 

0.07 


(0.47)5/4  (0.055)3/8  tanh  (( 


9  9  0.07 
0.47 


2.879  [m'1] 


Eq.  (10)  becomes 


iv  =  0.1352  1  + 


(i  +  lMD) 

V  9  (0.47)  / 


7  (0.07)  \  /  (1600)(2.935  107) 


0.47 


tanh((0.9589)(2.897)  d)\ 


2.897  d 

'1/2  fl/2  (  tanh 


dx  dy 


=  1503  x  4 


^2.778  yj x 


/o 


/o 


y/x2  +y2 


3.25 


dxdy  . 


This  integral  was  solved  with  the  Mathematica™  function  NIntegrate1  resulting  in  iv  —  17.37  [kN-s].  The  evaluation 
of  the  algebraic  Eq.  (11)  gives  iv  —  17.38  [kN-s],  a  difference  of  less  than  0.07%. 

1The  function  NIntegrate  is  documented  in  [6].  It  uses  the  Genz-Malik  adaptive  algorithm  to 
solve  the  multidimensional  integrals. 
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4.3  Square  Deflectors  Centered  above  the  Mine 

Consider  a  square  blast  deflector  of  width  l  centered  above  the  mine.  By 
substituting  x\  =  yi  —  1/2  and  x0  =  y0  =  — 1/2  into  Eqs  (13)  and  (14),  one  gets 


Si  = 


n 


«sC2 


erf2  (v^ti/C/2) 


(22) 


and 


S,  = 


«2 


(36  l2(2  +  6  k\  /4  C4)  e  K3i2<’"//2  +  96  7T  erf2  (ydyj  /  C/2) 


2  16/^C2  LV 

-/  C  (132  +  16  12  C2  +  4  l4  C4)  erf  ( vW  C/2)  e~^ * ^/4]  .  (23) 


4.3.1  Effect  of  the  Plate  Area 

When  the  area  A  =  /2  of  the  blast  deflector  increases,  the  total  vertical  impulse 
increases  also.  However,  if  the  standoff  distance  is  kept  constant,  the  impulse  will 
eventually  reach  a  limit  value  given  by 


lim  iv  =  Kq  1  +  ^  \  —  lim  (Si  +  S2) 

l->oo  \  9  Z  I  \  Z  l->  OO 


(1  +  ^)  vf  & 


7 r 


«3C2 


1  +  ^)  erf2  (yfcl  C/2) 


1.-3 

a3 


Ko  7T  \  / 1  6  k2 

I  I  1 


AsC2 


uc-* 

rvq 


Because  the  explosion  of  a  mine  under  a  deflector  of  infinite  area  is  the  worst 
possible  case,  this  equation  gives  the  upper  limit  of  the  impulse  on  any  blast  deflector 
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with  a  standoff  distance  of  z  meters: 


hnax 


/ 0.5857\  /  75\  fpE 

{—)  (1  +  9l)  VT 


(24) 


Example  2  Consider  the  parameters  of  Ex.  1  on  page  9  with,  this  time,  a  deflector  l  m  wide  by  l  m  long. 
Fig.  6  is  a  plot  of  Eq.  (11)  with  Si  and  S2  given  by  (22)  and  (23),  respectively.  The  maximum  impulse  computed 
with  Eq.  (24)  is  24.9  [kN-s]. 


FIGURE  6  -  Total  Impulse  as  a  Function  of  Plate  Width  for  Ex.  2 


4.3.2  Effect  of  the  Standoff  Distance 


The  total  vertical  impulse  iv  decreases  rapidly  with  the  standoff  distance  z. 


Example  3  Once  again,  consider  the  parameters  of  Ex.  1  on  page  9.  This  time,  the  standoff  distance  goes 
from  15  cm  to  4  m.  Fig.  7  shows  the  decrease  of  the  total  impulse  computed  with  Eq.  (11). 


5  OBLIQUE  BLAST  DEFLECTORS 


5.1  Specific  Impulse 


When  the  blast  deflector  is  not  horizontal  but  at  angle  0  <  a  <  n  from  the 
vertical,  as  in  Fig.  8,  the  resulting  specific  impulse  on  the  plate  at  location  P(x,y,z ) 
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FIGURE  7  -  Total  Impulse  as  a  Function  of  Standoff  Distance  for  Ex.  3 

is  also  a  function  of  a. 


FIGURE  8  -  Oblique  Blast  Deflector 


Let  ir  be  the  radial  specific  impulse  along  the  direction  of  the  target  point. 
From  Fig.  9,  the  specific  impulse  in  normal  to  the  blast  deflector  is  found  from  ir  by 


(in)  (clA)  =  (ir  cos  6)  (cos  0  clA ) 


or  in  =  ir  cos2  6.  Similarly,  iv  =  ir  cos2  /3.  By  combining  these  last  two  equations 
together,  one  gets  the  relation  between  the  known  component  of  the  specific  impulse 
iv  to  the  sought  component  in\ 


4  =  4  ^  •  (25) 

cos z  p 

Eq.  (25)  is  different  from  the  one  used  in  [7],  where  the  squares  of  the  cosines  are 
missing. 
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5.2  Total  Impulse 


For  a  blast  deflector  positioned  as  in  Fig.  10,  the  deflector  is  in  the  plane 
z  =  z\  +  —  x )  cot  a.  The  area  on  which  the  impulse  is  applied  is 

A  =  — — - —  (oblique  deflector  surface) 

sin  a 


and  the  element  of  integration  dA  is 


dA 


dx  dy 
sin  a 


The  total  normal  impulse  on  the  deflector  is  then 


rx  i 

rvi 

cos2  9 

)XQ  O 

'yo 

cos2  f3  sin  a 

in  =  0.1352 


dy  dx  (26) 
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where 


z  = 

A  —  x  cot  a 

(27) 

A  = 

Z\  +  Xi  cot  a 

(28) 

cos  9  = 

x  cos  a  +  z  sin  a 

(29) 

y/x2  +  y2  +  z2 

cos/3  = 

z 

(30) 

y/x2  +  y2  +  z2 

The  constant  A  is  the  standoff  distance  when  x  —  0.  By  combining  Eqs  (27)  to  (30), 
the  ratio  of  trigonometric  functions  in  (26)  becomes 


cos2  9  ( x  cos  ct  +  z  sin  a)2  A2  sin  a; 

cos 2  /3  sin  a  z 2  sin  a  (A  —  x  cot  a)2 


(31) 


FIGURE  10  -  Oblique  Blast  Deflector  over  the  xyz-Plane 


5.3  Approximate  Solution 


The  total  normal  impulse  in,  as  given  by  Eq.  (26),  cannot  be  analytically  inte¬ 
grated.  However,  with  the  approximation  (8),  an  algebraic  solution  can  be  obtained 
along  the  y  axis.  There  is  still  a  need  to  numerically  integrate  the  resulting  equation 
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along  the  x  axis  but  the  number  of  summations  drops  at  least  from  N2  to  N.  For  ex¬ 
ample,  if  a  grid  of  100  elements  (as  in  [7])  are  used  to  numerically  integrate  (26),  four 
elements  are  usually  enough  with  this  proposed  approximate  solution  (see  Ex.  5-7). 

5.3.1  Integration  Along  the  y  Axis 


16 


Example  4  Consider  a  blast  deflector  made  of  two  plates  of  size  0.5  m  by  1  m  each  and  attached  in  a  V-shape, 
as  shown  in  Fig.  11.  The  charge  is  a  TMA-3  anti-tank  mine  buried  with  3  cm  of  soil.  The  total  impulse  on  the 
deflector  is  equivalent  to  four  times  the  total  impulse  on  the  area  contained  in  the  region  x  E  [0,0.4]  and  y  E  [0,0.5]. 
Then,  from  Eq.  (35) 

/»0.4 

in  =  4  ko  I  p(x)  dx 

Jo 

and  from  Fig.  11,  yo  =  0,  yi  =  0.5,  cotct  =  —0.75  (or  a  =  126.9  deg)  and  2:  =  0.47  +  0.75#.  From  Ex.  1  on  page  9, 

5  =  0.07  [m] 

p  =  1600  [kg/m3] 

E  =  2.935  107  [J] 

A  =  0.055  [m2] 


and 


c  = 


0.2077 


(0.47  +  0.75 1)5/4  tank 


0.0604 


v  3/2 


^  0.47+0.75  x  7 

The  function  f  =  <^(z(x))  for  0  +  .T  +  0.4  is  plotted  in  Fig.  12  and  is  well  approximated  by  the  line 


The  graph  of  p(x)  is  shown  in  Fig.  13. 


C  ~  2.887 +  0.872  a;  . 


FIGURE  11  -  V-Shaped  Blast  Deflector  Centered  Above  the  Mine 


5.3.2  Integration  Along  the  x  Axis 


Let 


£ 


2x  —  xq  —  x\ 

J 


Xi  -  x0 
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FIGURE  12  -  Graph  of  C(z(z))  for  Ex.  4 


FIGURE  13  -  Graph  of  Eq.  (36)  with  the  Parameters  of  Ex.  4 
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then  Eq.  (35)  becomes 

i„  =  Kon_^£  |’p(lMwa)  dg  .  (37) 

This  equation  can  be  numerically  integrated  using  a  Gaussian  quadrature.  A  4-point 
Gauss-Legendre  integration  is  usually  enough  to  integrate  a  smooth  function  such  as 
p(x)  in  Ex.  4  on  page  16  (Fig.  13). 

The  total  normal  impulse  on  an  oblique  blast  deflector  can  then  be  approxi¬ 
mated  by 

(38) 

where  G;  and  Wi  are,  respectively,  the  abscissas  and  weights  of  the  JV-point  Gauss- 
Legendre  quadrature  formula.  Values  from  [8]  of  G  and  wt  for  N  =  4  and  6  are  shown 
in  Table  I. 

TABLE  I 

Gauss-Legendre  Abscissas  and  Weights 


N 

i 

G 

Wi 

4 

1,2 

±0.33998  10435  84856 

0.65214  51548  62546 

3,4 

±0.86113  63115  94053 

0.34785  48451  37454 

6 

1,2 

±0.23861  91860  83197 

0.46791  39345  72691 

3,4 

±0.66120  93864  66265 

0.36076  15730  48139 

5,6 

±0.93246  95142  03152 

0.17132  44923  79170 

Example  5  The  total  impulse  on  the  V-shaped  deflector  of  Ex.  4  on  page  16  is  given  by  Eq.  (38)  as 
in  —  4  x  0.0236  [W!  (p(0.1320)  +  p(0.2680))  +  w3  (p( 0.0280)  +  p(0.3722))] 

=  6.97  [kN-s]  . 

TM 

The  integral  (26)  was  solved  with  the  Mathematica  function  NIntegrate  resulting  in  in  —  6.99  [kN-s],  a  difference 
of  0.3%. 

5.4  Deflectors  Centered  on  the  x  Axis 

Let  l  be  the  length  of  the  deflector  along  y.  Because  the  deflector  is  centered 
on  the  x  axis,  the  integral  along  y  for  — 1/2  ^  y  ^  1/2  is  twice  the  value  of  the  integral 
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on  the  interval  0  ^  y  ^  1/2.  Then  Eqs  (33)  and  (34),  used  in  the  calculation  of  p(x), 


can  be  written  as 


T^x^)  =  erf  (v^aC l/2)  \  — tv  e 


—  K3  C2  X 2 


3-k3  C2  X2 


r2 (x,z)  =  2  '  7/2  1  [-«3/2  c5  ^5/4  -  «3/2  c3  (5  +  6K3C2a;2)  /2 

8«s  C  ( 

—3  (l  (^5  +  6  k3  £ 2x 2  +  4  (4  a;4)  e~K3  ^ ,2//4  + 

y/rr  ^15  +  18  k3  £ 2  x2  +  12  k2  ^4a;4  +  8  ^  C6^6)  erf(y^£  Z/2)  \  .  (40) 


6  VERTICAL  BLAST  DEFLECTORS 


6.1  Specific  Impulse 


When  the  blast  deflector  is  vertical,  a  =  0  or  n  and  the  deflector  is  in  the 
plane  x  =  x0  (Fig.  14).  The  deflector  is  in  the  region  y  e  [yo,y  1],  7  G  [zq,  Zi]  and  an 
approach  similar  to  Chapter  5  is  applied.  The  ratio  of  cosines  in  Eq.  (25)  becomes 

cos2d  (Xq\2 
cos2  (3  \  z  ) 

and  the  horizontal  specific  impulse  on  the  vertical  blast  deflector  is 

ih  =  iv  (^7)  •  (41) 


6.2  Total  Impulse 


The  total  impulse  on  a  vertical  blast  deflector  is  given  by 

■" = °'i352fr  1(t)2  f1 + ¥z) n(d)  0\ dydz 
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,  z 

n 


FIGURE  14  -  Vertical  Blast  Deflector 


where  d  = 


6.3  Approximate  Solution 


By  applying  the  method  of  Section  5.3,  the  total  impulse  on  a  vertical  blast 
deflector  can  be  written  as 


where 


(43) 


def  X20  v fpE 
«(*)  =  ^5/2 


ri(a;o,z)  +  r2(iCo  ,z) 


(44) 


and  ri(a;o,*),  r2(a;o,^)  are  given  by  Eqs  (33)  and  (34)  with  x  =  xQ.  Eq.  (43)  can  be 
approximated  by 


(45) 


where  £*  and  Wi  are  given  in  Table  I. 
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Example  6  Consider  a  square  vertical  blast  deflector  of  area  A  =  1  m2  and  centered  on  the  x  axis.  Let  us 
predict  the  total  impulse  when  xo  E  [0,  oo)  and  the  remaining  mine  and  plate  geometries  are  as  those  of  Ex.  1  on 
page  9.  Eq.  (44),  with  the  use  of  Eqs  (39)  and  (40)  for  a  symmetric  problem,  is  plotted  in  Fig.  15.  This  function 
was  then  integrated  with  a  4-point  quadrature  (Eq.  (45))  and  plotted  in  dashed  lines  in  Fig.  16;  the  solid  line  is  the 

TM 

solution  of  Eq.  (42)  directly  from  Mathematica  .  Discrepancies  between  the  two  curves  arise  from  the  difference 
between  f((,d)  and  g(C,d)  for  C,d>  2.5  (Fig.  3). 


FIGURE  15  -  Graph  of  Eq.  (44)  for  Ex.  6 


FIGURE  16  -  Total  Impulse  as  a  Function  of  xq  for  Ex.  6 


Ex.  6  above  illustrates  how  the  specific  impulse  modeled  by  Westin  et  al.  in  [4] 
is  sensitive  to  the  lateral  distance  d  =  \Jx2  +  y2  from  the  center  of  the  mine.  In  Ex.  3 
on  page  11,  the  effect  of  the  impulse  is  still  being  felt  4  m  above  the  mine  (Fig.  7). 
However,  in  Ex.  6,  the  maximum  impulse  is  at  a  lateral  distance  of  50  cm,  and  drops 
rapidly  after  that  (Fig.  16).  In  both  cases,  the  plate  was  1  m2.  That  is  a  consequence 
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of  the  assumptions  that  the  loading  produced  from  a  landmine  explosion  is  mainly 
the  result  of  detonation  products  and  soil  ejecta,  which  are  propulsed  upwards  like  in 
a  cone,  and  that  blast  overpressure  plays  only  a  secondary  role. 


Example  7  Consider  the  V-shaped  deflector  in  Fig.  17  where  7  varies  from  0  (a  vertical  plate)  to  180  deg  (a 
horizontal  plate).  The  two  1  m  long  plates  are  centered  above  the  mine  and  the  total  deflector  has  a  presented  area 
A  =  1  m2  on  the  xy  plane  for  7  >  0.  The  total  impulse  on  the  deflector  is  twice  the  impulse  on  one  single  plate. 
Then, 


7  G  [0, 180]  [deg] 

a  =  7T  ( 1 - —  )  [rad] 

V  360/  L  J 

£o  =  0 

x\  =  0.5  [m] 

zo  =  A  =  0.47  [m] 

z\  =  X  —  0.5  cot  a  [m] 

V  £  [—0.5,  0.5]  [m] 

and  the  impulse  is  twice  the  value  predicted  by  Eqs  (11),  (38)  or  (45).  The  total  impulse  for  selected  angles  7  are 
given  in  Table  II  and  plotted  at  Fig.  18  along  with  Eq.  (38),  using  N  =  4. 


FIGURE  17  -  V-Shaped  deflector  of  Variable  Angle  7 
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TABLE  II 

Impulse  on  the  Variable  V-Shaped  Deflector  for  Ex.  7 


7  [deg] 

a  [rad] 

Equation  used 

Impulse  [kN-s] 

0 

71 

(45) 

0 

22.5 

15  tt/16 

(38) 

0.5 

45 

7  7t/8 

(38) 

1.7 

90 

3  7t/4 

(38) 

5.5 

135 

5  7t/8 

(38) 

10.6 

180 

tt/2 

(11) 

17.4 

FIGURE  18  -  Impulse  as  a  Function  of  the  Angle  7  in  Fig.  17 
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7  CONCLUSIONS 


Algebraic  equations  were  derived  to  predict  the  total  impulse  on  a  flat  blast 
deflector  at  an  arbitrary  orientation.  Furthermore,  a  very  simple  equation  predicting 
the  maximum  impulse  at  a  given  standoff  distance  was  also  derived.  These  equations 
can  easily  be  coded  in  either  a  spreadsheet  or  a  computer  program. 

As  they  stand,  these  equations  are  an  aid  to  engineers  who  need  to  estimate 
the  loadings  to  be  expected  for  the  design  of  blast  deflectors.  However,  more  research 
is  required  to  model  and  predict  the  effects  of  the  impulse  produced  from  a  mine 
explosion  under  a  wheel  combined  with  blast  deflectors  loadings. 
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APPENDIX 

The  Error  Function 


The  function  erf  (re)  is  the  error  function  defined  by 


erf  (re) 


A  =  Jp(A2)  if  x  ^  0 

VWo  (-P(4,re2)  if  re  <  0 


where  P{\,  re2)  is  the  incomplete  gamma  function.  Reference  [9]  proposes  a  Chebyshev 
fitting  for  erf  (re)  with  a  fractional  error  everywhere  less  than  1.2  10  ' : 


erf  (re)  ~ 


re  ^  0 
re  <  0 


with 


9  =  t  exp  (—re2  -  1.26551223  +  t  (1.00002368  +  t  (0.37409196  + 

t  (0.09678418  +  t  (-0.18628806  +  t  (0.27886807  +  t  (-1.13520398  + 

t  (1.48851587  +  t  (-0.82215223  +  0.17087277 1)))))))) 


and 


